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ABSTRACT 



An evaluation is performed of a radiation model for the 
Naval Postgraduate School primitive-equation numerical wea- 
ther prediction system. The model employs empirical expres- 
sions for atmospheric absorptivity, scattering-reflectivity, 
cloud-reflectivity and earth-surface reflectivity to compute 
solar insolation absorbed at earth and in the key atmos- 
pheric layers. The terrestrial cooling effect at earth and 
in these same key atmospheric layers is formulated using re- 
cent empirically-derived emissivities for the effects of 
both water vapor and 

Mean seasonal atmospheric soundings for the Northern 
Hemisphere are utilized for testing the model. In addition, 
application of atmospheric boundary-layer modeling permits 
determination of the surface-layer turbulent transports (in 
the vertical) of sensible and latent heat at the earth’s 
surface. To evaluate the validity of the radiation model, 
heat budgets are compiled for the earth-atmosphere system, 
the atmospheric column and the earth’s surface. 
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INTRODUCTION 



I 

f 

I . 

The present study was to evaluate a radiation package 
for the Naval Postgraduate School primitive equation pre- 
diction model. The radiation package is primarily an out- 
growth of the model used in the UCLA general circulation 
model (Gates et al, 1971) with modifications suggested 
from the NCAR model (Oliger et al. 1970) and from the Geo- 
physical Fluid Dynamical Laboratory (Manabe and Strickler , 
1964) . 

The model consists of five a-levels ( cr = P/tt where TT 
is the surface pressure) for the purposes of data input 
as shown in Figure 1. Cloud cover was treated by the use 
of a single cloud layer extending from O = 0.8 to 0.4. 

In the computations for both pressure-scaled water-vapor 
and CO 2 absorber masses, a more realistic method was in- 
troduced to extrapolate accurate water-vapor values above 
levels normally available from radiosonde' data (P ~ 300 mb) . 
A system for the integrations of net terrestrial flux was 
introduced by using empirically derived absorpt ivit ies 
(Sasamori, 1968) based on the radiation chart of Yamamoto 
(1952), for water-vapor and 00^* Evaporation and sensible 
heat fluxes were derived from a formulation for surface 
layer heat transfer by Langlois and Kwok (1969). 

Northern hemisphere climatological atmospheric data 
from London (1957) were utilized to evaluate the model. 
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Pressure, temperature, and relative humidity data were 
divided into the four seasons and available over ten de- 
gree latitude bands from the equator poleward. 

The development of the terrestrial fluxes is given 
in Section III and the solar radiation model is given in 
Section IV. The results of the computed long and short 
wave fluxes were then compared to those of Budyko (1955), 
London (1957) and Katayama (1967). 

Finally, to evaluate the adaptability of the radiative 
package to a general circulation model, three heat bud- 
gets were attempted. The budgets were for the earth-at- 
mosphere system, the atmospheric column and the surface. 
Since heat transfer calculations due to evaporation and 
sensible heat are more easily made over ocean surfaces, 
London's atmospheric data was fitted to climatological 
ocean surface temperature and wind data from the United 
States Weather Bureau’s Atlas of Climatic Charts of the 
Oceans (1938). Surface layer evaporation and sensible 
heat flux computations were not conducted poleward of 60N 
because of the possibility of ice cover. Rather, com- 
bined evaporation and sensible heat rates were taken from 
Vowinckel and Taylor (1964). 
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k-Level 



0-Level 



Layers 

T 

( 0 , 2 ) 

i 

(2,4) 

I 

(4,6) 

\ 

( 6 , 8 ) 

i 

( 8 , 10 ) 

JL 



0 0.0 

1 Au , Ac , q , T 0.1 

2 T 0.2 

3 Au , Ac , q 0 . 3 

4 T 0.4 

5 Au , Ac , q 0 . 5 

6 T 0.6 

7 Au,Ac,q 0.7 

8 T 0 . 8 

9 Au , Ac , q 0 . 9 

10 T 1.0 



Figure 1. Representation of the five-layer sounding 
model used for radiative transfer calculations-. Tem- 
perature, mixing ratio, pressure-scaled water-vapor 
(Au) and CO 2 (Ac) absorber mass increments are shown 
qt the appropriate O-levels and k-levels (k = lOo) . 
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II. DATA PREPARATION 



A. TEMPERATURE AND MOISTURE TERMS 

Northern hemisphere vertical and latitudinal distribu- 
tions of pressure, temperature, and relative humidity were 
taken from the climatological soundings of London (1957). 
The data was compiled in four mean seasonal tables. The 
latitudinal variation consisted of 10 degree latitude bands 
measured from the equator northward. In the vertical, data 
was distributed at 1 km intervals from the surface to the 
base of the stratosphere. 

With actual radiosonde data, continuous information 
would be available above the tropopause to at least the 
0.1 sigma level. The 0 = 0.1 level was the highest level 
at which the water-vapor mixing ratio was considered to be 
recoverable from the fields of q. London’s radiosonde data 
was continuous only to the tropopause. Since the tropo- 
pause varies in height lat it udinally as well as seasonally, 
adjustment of the vertical distribution of T and q had to 
be made to insure that the data would be available to O = 
0.1 level. The simplest solution for the stratosphere was 
to assume isothermal conditions with constant relative 
humidity to the O = 0.1 level for any vertical sounding 
that ended below 0 = 0.1. 

It was necessary to interpolate from the distribution 
at successive km levels in the vertical (as listed in 
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London, 1957) to the specific sigma levels of Figure 1. 
Temperatures were carried at the 1.0, 0.8, 0.6, 0.4, 0.2 
and 0.1 G-levels. Since surface data was available in 
London’s soundings, was set equal to the listed sur- 

face temperature rather than requiring an extrapolation 
of the type 

T 1q = 0.5 (3T 9 - T ? ) (1) 

Temperatures at other sigma surfaces, say at level k, were 
interpolated using a three-point Lagrangian formula of the 
form 

TOO - (T 2 )(P k -P 1 )(P k -P 0 )/(P 2 -P 1 HP 2 -P 0 ) 

+ (T 1 )(P k -P 2 )(P k -P 0 )/(P 1 -P 2 )(P 1 -P 0 ) 

+ (T 0 )(P k -P 1 )(P k -P 2 )/(P 0 -P 2 )(P 0 -P 1 ) (2) 

where it is assumed that the data level (Pq,Tq) lies above 
the G-level in question and (^2*^2^ ^-^ es below that level. 

Mixing ratio values were calculated with London’s data 
us ing 

q ** RH (0.622e ) /P (3) 

s 

The q-values were then interpolated to G^-levels using the 

Lagrangian formulation of (2) applied to q’s at known 

levels. The q-values were centered at the odd G-levels 
k 

of 0.9, 0.7, 0.5 and 0.3. These values are then considered 
to be the layer means between the even G-levels. 

Accurate humidity information is in general not avail- 
able above 0 = 0.3, especially in the stratosphere. To 
find a mean q-value representative of the layer (0,2) a 
different formulation was used. 
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Smith (1966) has shown a strong correlation between 
q(P) and the surface value through the extrapolative 

formula 



q(P) = f P q 

q io ' p io 



(4) 



London’s soundings were tested to obtain a best-fit be- 
tween In q and In P of the form 



q(P) = 

q 3 



(^-) X 



(5) 



which is similar in form to that suggested by Smith, The 
values of q(P) resulting from (5) for O = 0.9, 0.7, 0.5 
and 0.3 were used to solve for the profile parameter X. 
This parameter varied with latitude and season as shown 
in Table I. The relationship (5), with appropriate 
X-values, was then used to extrapolate q at a = 0.1. The 
resulting q^-values were then considered the mean 
q-values for each layer (0,2). 

The procedure just described to derive q^ was re- 
garded as justifiable in view of the following reasons: 
-.(i) The residuals [q(est) - q(obs)]^ were small for 
all k >. 3 , and 

(ii) The extrapolated q^ corresponds to the proper 
order of magnitude as listed in model atmospheres (Valley, 
1965) at pressures P - 100 mb. 
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TABLE I. Seasonal and latitudinal values of the profile 
parameter X of Equation (5). 



Latitude 


Winter 


Spring 


Summer 


Fall 


0-10 


3.1438 


3.0178 


3.2885 


3.0202 


10-20 


3.1182 


3.3024 


3.4381 


3.2555 


20-30 


3.1205 


3.3978 


3.4272 


3.2374 


30-40 


2.6437 


3.3838 


3.3563 


3 .0873 


40-50 


3.0002 


3 . 5178 


3 . 2480 


3 . 1259 


50-60 


2.7210 


3 . 2382 


3 . 2031 


3 .1846 


60-70 


2 . 3267 


2.7478 


3.0684 


2.8925 


70-80 


1.9671 


2.1953 


2.7886 


2.4617 


80-90 


1.1743 


1.7831 


2 #8434 


2.1262 
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B. CLOUD MODELING 



In both the solar and terrestrial radiative transfer 
processes, it is essential to know the fractional cloud 
cover (CL) at each grid point. Clouds act as almost perfect 
blackbody radiators for infrared radiation, and as efficent 
reflectors for solar radiation. 

There is a difficulty in parameterizing different 
cloud types and heights. In this model, as in the UCLA 
general circulation model, only a single cloud layer was 
parameterized for the purposes of inclusion in the vertical 
radiative transfer model at gridpoints. As is done in 
most other numerical prediction schemes involving cloud 
amount, CL, use is made of Smagor insky 1 s (1960) algorithm 
(applicable at approximately 700 mb) 

CL = 2.0 (e/e e ) - 0.7 (6) 

where 0 <, CL 1.0. When a cloud amount CL > 0 existed, 
it was assumed to extend between O = 0.8 to O =0.4. 

v 

Comparison of the calculated values of CL computed 
from (6) and those tabulated from seasonal climatological 
estimates by Telegadas and London (1954) and Seinde (1954) 
showed a marked difference in the mid-latitudes as seen 
by comparison of Tables II and III. Since other climato- 
logical estimates of cloud cover, namely Sellers (1965), 
compared best with the data tabulated in Table II, climato- 
logy was used in this model rather than the results of 
Table III. 
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TABLE lx. The seasonal distribution of total cloudiness 
(percent cloud cover) tabulated by London (1957). 



Lat N. 


Winter 


0-10 


47 


10-20 


36 


20-30 


38 


30-40 


50 


40-50 


59 


50-60 


63 


60-70 


58 


70-80 


47 


80-90 


40 



Summer 


Fall 


54 


53 


49 


48 


42 


41 


41 


46 


55 


56 


63 


66 


66 


70 


69 


70 


64 


60 



Spring 

51 

42 

42 

52 

59 

62 

60 

59 

55 
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TABLE III. The seasonal distribution of total cloudiness 
(percent of sky cover) calculated by the Smagor insky 
algorithm . 



Lat N 


Winter 


Spring 


Summer 


Fall 


0-10 


30 


40 


56 


46 


10-20 


12 


20 


. 48 


36 


20-30 


10 


10 


30 


20 


30-40 


18 


14 


22 


16 


40-50 


34 


26 


30 


26 


50-60 


42 


34 


42 


42 


60-70 


46 


40 


50 


46 


70-80 


46 


38 


50 


48 


80-90 


34 


28 


50 


40 



17 



Even though London’s cloud-cover climatology was 
accepted as valid here, it should be realized that London’s 
results included overlapping clouds at multiple levels, 
which was impossible with the parameterization using (6). 

An improved cloud model permitting cloud existence at 
two or more levels in the vertical should be attempted 
in future radiative heating systems. At the present time, 
the added complexity of the radiative physics is considered 
too costly in computer time for inclusion in the present 
model. However, it should be noted that Arakawa et al. 
(1972) have already developed the concept of a three-layer 
radiative cloud model for a 1972 version of the UCLA 
general circulation model. 

Another reason for the choice of the climatological 
CL must be considered. With the use of climatology, 
transient synoptic scale perturbations and disturbances 
in the soundings are smoothed or eliminated. Therefore, 
the day-to-day variations of relative humidity tend to be 
smoothed out during the season. The resulting smoothed 
relative humidity distributions thus tended to be less 
than their normal peak-values and give too small a magni- 
tude of CL when employing the Smagorinsky algorithm, 

C. ABSORBER-MASS DEFINITIONS: PRESSURE SCALED 

Since mixing ratio values have now been computed at 
the five odd a-levels (Figure 1), it is possible to deter- 
mine the pressure-scaled water-vapor absorber masses. 



18 



Au^> associated with the five corresponding layers. This 
is done in the model in such a way that the additive 
s cheme 



u ( 6 , 10) = u ( 8 , 10 ) + u ( 8 , 6 ) (7) 

holds. The layer water vapor absorber mass, u(2k, 2k-2) , 

has been centered at odd levels (2k-l) . The scaled 
absorber mass is calculated using 



Au 



q 2k-l AP / P 2k-1,0. 



(■ 



-) 



72 



( 8 ) 



2k-l g ‘o 

which, as has been noted, is the absorber mass in the 

layer (2k, 2k-2) . AP is a constant in the analysis scheme 

equal to 0 . 2 tt . The mixing ratio value, ^k-l* t * le 

mean centered value in the layer (2k, 2k-2) . The pressure 

0.72 

ratio factor ( P 2k-l^ P 0^ * the sca ^-^ n S factor for 

collisonal line-broadening after Moller and Raschke (1964). 
This latter formulation disposes of the necessity of 
temperature scaling in the usual Lorentz line-broadening 
formula (c.f., Danard, 1969), where the pressure-ratio 
scaling is linear and a temperature-ratio scaling factor 
is indicated by (T/Tq) times the N.T.P. mass within 

a layer. The subscript "0" denotes an N.T.P. condition 
where Pq = 1013.25 mb, T n = 273.16°K, and p n is the N.T.P. 



value of density. 



0 



0 



If (8) is integrated from = (0*2k) Pq to P = 7T , 

the result is 

17 2 P 0 n ~ _ 17 2 17 2 

U 2k = (7r/P 0 ) ( 1. 7 2g ^ Z q 2n-l ^°2n " a 2n-2 ] 

n= j 
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where by definition 



, _ P 2k P 0 

°k = P 2k /7T = ~ 



( 10 ) 



For Roller and Raschke (1964) proposed pressure- 

scaling similar to that of water-vapor but with a pressure- 
ratio exponent of 0.65. For simplicity, a single 
pressure-scaling exponent of 0.72 was utilized for both 

cases. The form of the layer-reduced absorber mass of CO^ 

2 

in N.T.P. cm/cm in the layer (2k, 2k-2) is 

Ac 2k _i = 3.14 x 10~ 4 ^AP/gp 2k-l ) 0.72 (11) 

The process of integration from to 

formulation 



'2k 



3.14 x 10 
1.72 



-4 P 1.72 

— H[1 - ] 

0 



( 12 ) 



where H = 7.995 x 10 cm is the height of the homogeneous 
atmosphere in N.T.P. cms . Adjusting (12) to a O-coordinate 
system with P = u at the surface yields 



- 2k 



3.14 x 10 
1.72 



-4 



H ( f o ) 1 - 72 cx- o 2k 



1.72 



] 



(13) 



Table IVshows the seasonal distribution of the total 
pressure-scaled water-vapor absorber masses, u(0,10). As 
expected, the summer season contains the largest absorber 
mass with the least in winter. Spring and fall values are 
intermediate between the maximum and minimum values. Com- 
putations for the CO 2 absorber-masses are dependent upon 
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TABLE IV. Seasonal distributions of total pressure-scaled 



water-vapor 


masses , u (0 


,10) , in gm 


“2 

cm^ . 




Latitude 


Spring 


Summer 


Fall 


Winter 


0-10 


3.878 


4.287 


4.023 


3.624 


10-20 


3.271 


4.316 


‘ 3.826 


2.749 


20-30 


2.417 


3 . 601 


3.055 


1.993 


30-40 


1.607 


2.897 


2.133 


1.279 


40-50 


1.099 


2 .292 


1.426 


0 . 747 


50-60 


0.703 


1.793 


0.932 


0.449 


60-70 


0.437 


1.428 


0.598 


0.256 


70-80 


0 .289 


1.045 


0 .471 


0.163 


80-90 


0.192 


0 .901 


0.319 


0.105 
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the total mass of air in a unit column, and hence on the 

surface pressure. Since the variation in surface pressure 

averaged across all latitude bands was less than 10 mb, 

the CO^ absorber-masses were nearly constant [c(0,10) = 

2 

125 cm/ cm ] . 
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III. TERRESTRIAL RADIATIVE TRANSFER 



A. THEORETICAL AND EMPIRICAL BACKGROUND 

Virtually all of the radiation packages in present use 
with large scale circulation models in the United States 
[cf , , Manabe and Strickler (1964), Oliger et al. (1970), 
Arakawa et al . (1969)] make use of Yamamoto's (1952) water- 

vapor flux emissivities as a basis for their long-wave 
radiative flux calculations. The water-vapor flux emissi- 
vity is normally thought of as a function of both water- 
vapor scaled mass and temperature. Yet, in the normal range 
of temperatures (T >L 220K), the emissivity is virtually 
independent of temperature. At temperatures T < 220K, the 
wave-averaged flux emissivity is defined after Yamamoto 
(1952) by 

00 dB 

1 11 " t FV <1 V0 U)I 7 t-< T - 2 2° K > ** 

e(u.I) - - 2 S <1 4 > 

. dB 

O ir <V>T) dv 

Here the numerator has its temperature dependent terms 

3 

fixed at 220K, while the denominator equals 4St(T ) even 
when T<220K. The temperature affects the population of 
different rotational energy levels within the so-called 
"rotational 11 band. At low temperatures, Planck's function 
shifts maximal radiant energy to the longer wavelengths 
in the rotational band, according to Wien's displacement 
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law o f 



X = 2897 /T (15) 

max ' 

where X is in microns. Thus with T < 220K, the highly 
absorbed rotational lines are the only contributors to 
water-vapor absorptivity. This leads to an increase of 
absorption near the center of the rotational band and thus 
to an increase in mean absorptivity with decreasing tem- 
perature. This phenomenon is depicted graphically in 
Figure 2 after Sasamori (1968). 

Yamamoto’s radiation chart has been the basis of most 
recent radiation calculations. Sasamori (1968, 1970) used 

Yamamoto’s radiation theory to write empirical formulas 
for all of the emissivities needed for long-wave radiative 
flux calculations. These empirical formulas have been 
shown by Sasamori to yield adequate accuracy and were 
adapted for this numerical study. The temperature indepen- 
dent emissivity for water-vapor, from Sasamori, is given by 
e (u) = 0.24 log 10 (u + 0.010) + 0.622 (16) 

This formulation was used for T >_ 220K. For T < 220K, 
Sasamori gives the temperature-dependent part of the emis- 
sivity function by water vapor as 

€(u,T)-(8.34T 0 - 3531os 10 u - 0 - 44 )u-°- 34551o 8 10 u '- 705 (17) 

When absorption by CC^ and water-vapor are considered 
jointly, the product law of transmissivities must be 
employed. In doing this, the transmissivity of the two 
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e (u) 




u (gm 



cm ^ ) 



Figure 2. Temperature dependent part of the emissivity 
with varied amounts of scaled water-vapor. 
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absorbers acting jointly may be expressed as 



T p (u,c) = T p (u)T F (c) 

The flux integral as formulated by Yamamoto is 



( 18 ) 



F 



u(T) 

/ e(u,T) 



u=0 , T 



Re f 



dB 



(19) 



where u(T) is the water-vapor sounding mass measured away 
from the reference level. Equation (19) then becomes with 
the inclusion of CO^ 

u , c (T) 

F = / [l-T_(u)T_(c) ]dB (20) 

u=0,Q=0 F F 

or 



u >c 

F = / [e ( u , T) + Ae] dB (21) 

u = 0 , c = 0 

The Ae term represents the effective emission due to CO^ 
not already counted in the water-vapor flux calculations. 
The added emission is expressed on the Yamamoto chart by 
an increase of emissivity at all points of the radiative 
sounding (except at the reference level points). The Ae 
term includes only the 15ym band emitted flux of CO^ which 
has not already been counted in emission by overlap with 
the rotational band of water vapor. Sasamori considered 
this CO^ band emissivity to be temperature independent 
in developing the empirical formulation 

Ae = 0.07262[l-0. 62556 (u+0 .286) 0 * 26 ] [log c +1.064] 



( 22 ) 
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In the following flux computations, the emissivity, e, 
will include both the effects of CO^ given by (22) and that 
of water vapor given in (16) and (17). 



B. TECHNIQUES OF NET FLUX COMPUTATIONS 

Schematic representations of the Yamamoto radiation 
chart were utilized for considering the integration of the 
net long wave flux calculations. In Figure 3, the hatched 
area is a representation of downward flux at the surface. 
The area may be computed using the trapezoidal summation 
rule as follows: 

F+(10) = 0.5{ [e(8-10) (B 10 -B 8 ) ]+[e(8-10)+e(6-10) ] 
(B g -B 6 )+[e (6-10 )+e (4-10) ] (B g -B 4 )+[e(4-10) 

+ e(2-10) ] (B 4 -B 2 )+[e(l-10)+ e(2-10)] 

B(T ) 

(B 2 -B 1 )+ 2 / X e [ (0-10) ,T]db} 



4 

Here B - St * T is the Stef an-Bolt zmann blackbody flux 
arid e(8-10), for example, represents the emissivity com- 
puted from Sasamori’s empirical equations using the water 
vapor and CO^ present in the layer (8-10). The term 
S = e(u,T) + Ae is solved by e (u,T) given by(17)and Ae by 
(22). The subscripts on the S t e f an-Bol t zmann fluxes (i.e., 
B 10~ B 8^ are i^antif ers f° r t ^ le layer contribution from 
which the flux-emission takes place. 
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Figure 3. Schematic area-depiction of net flux at 
level 10. 
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e dB , 



B(T ) 

The final integral of (23), / 

B = 0 



may be evaluated by the integrated result 



St * T 4 { £ (u , T. ) (-- — , _ 

1 o 1 0.353 log-,, u + 3.56 

& 10 o 



) + 



0.07262 [1-0.62556 (u q + 0 . 286 ) 0 * 2 6 ] [ 1 og 1() c o + 1.064]} 



(24) 



In (24), it is understood that the integration spans 
the temperature range (0,T^) along the final isopleths 
u q and c q (e.g., Figure 3). The net flux at level 10 is 
given by 

F io - B io - F *< 10 > <25 > 



where the asterisk superscript denotes F(up) - F(down) at 
the level . 

The net flux at level 6 is represented schematically 
by B 1q minus the hatched area in Figure 4, The numerical 
modeling follows (23) where 

F 6 * = B 10 - 0.5 {e (6-8) (B g - b 6 ) + [ e (6-8)+ e(6-10)] 
(B 10 - Bg) + e (4-6) (Bg- B 4 ) + [£(4-6)+ £(2-6)] 
(B 4 - B 2 ) + [£(2-6) + e( 1-6)](B 2 - B x ) + 

B(T ) 

2 / e [ (0-6) , T] dB} (26) 

B = 0 
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Figure 4. Schematic area-depiction of net flux at 
level 6. ' 
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Figure 5 is the schematic for net flux at level 2. 
Following the previous computations 

F 2 * = B 1Q - 0.5 {e(2~4) (B 4 -B 2 )+[e(2-4)+ e(2-6)] 

( B 6 - B 4 ) + [e(2-6)+ e(2-8)](B g - B g ) + 

[ e (2-8) + e(2-10) ] (B 1q - B g ) + e(l-2)(B 2 - B ± ) 

B (T ) 

+2 / A e [ (0-2) , T] dB} (27) 

B = 0 

& * * 

The foregoing net fluxes F_ , F, , F n are clear sky 
calculations only, A modification due to cloudy skies is 
made. The c lo ud - co v er ed fraction is given by CL and the 
clear sky fraction is (1-CL) . The cloud-covered fraction 
is considered an overcast which extends vertically from 
level 8 to level 4. Figure 6 shows a schematic flux dia- 
gram, after Yamamoto (1952), of net fluxes at levels 10 
and 2, when only the overcast sky fraction is considered. 
The horizontally hatched area is now the net flux at level 
10. The net flux at level 6 is zero since a layer of ap- 
proximately 50 meters in depth with average cloud droplet 
content may be regarded as a blackbody (Brunt, 1939) at 
the mean temperature of the cloud layer. The net flux at 
level 2 is the slant hatched area in the diagram. 

The net flux at level 10, with overcast skies is given 
by 

F 10* = (B 10 “ V (1 " 0 . 5e (8-10) ) (28) 

Normally (28) represents a cooling effect at the earth's 
surface. At the same time, the net flux at level 8, or the 
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Figure 5. Schematic area-depiction of net flux at 
level 2. 
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Figure 6. Schematic area-depiction of net flux at 
reference levels 10 and 2 when an overcast cloud 
extends between levels 8 and 4. 
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cloud base, is an equal magnitude (Fig. 6) so that the 
layer (8-10) has zero flux divergence. However, the cloud 

* 

base receives a net flux warming equal in magnitude to 

The net flux at level 2 for the overcast case is given 
by 

F 2 * - B 4 - 0.5{e(2-4) (B 4 - B 2 )+ e(l-2)(B 2 - B^ 

B(T ) 

+ 2 f El (0-2) , T] dB} 

, B=0 ( 29 ) 

* 

This net flux, F^ , is much smaller than that for clear 
skies because of the lead term on the right side of 
(29) as contrasted with B^q in (27). 

The net fluxes, computed at the three reference levels, 
must now be clo ud -we igh t e d to yield final results for net 
long-wave flux transfer across the interface by 

F 10* = I’?0 (28) ] CL + lF 10 (25)] (1-CL) (30) 

F 6 * = [F*(26)] (1 - CL) - (31) 

F 2 * = [F*(29)] CL + £ F *( 2 7 ) ] (1-CL) (32) 

The temperature-dependent flux integral denoted by 

/ £ ( U Q > T^) was evaluated for an entire range 

of final T-. and u values which 
1 o 

could occur in the atmosphere. These results are listed 
in Appendix Table A-l. 

C. DIVERGENCE OF NET FLUX AND RESULTANT COOLING RATES 

The resultant, cloud weighted, seasonal distributions 
of net flux at the three reference levels of a = 0.2, 0.6, 
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and 1.0 are listed in Tables V(a-d). The 24-hour tempera- 
ture changes between the reference levels due to terrestrial 



radiation alone are also listed. The mean radiative tem- 
perature change centered at the midpoint of the column was 
computed by the first law of thermodynamics so that the cool- 



ing rate for the layer (6,10) is 



H _ * * 1000 

3t (F 6 F 10 )/0,47T g c p 



and for the layer (2,6) is 



3T rv * v *\ /n i-n 1000 

3T ( F 2 - F 6 )/0 * 47T ~ T % ’ 



(33) 



(34) 



where 7T is surface pressure in mb. 

To compare the accuracy of the long-wave calculations, 
the meridionally averaged mean values for the fluxes at the 
levels a = 0.2 and 1.0 were computed. The result was repre- 
sented by the percentage of the meridionally averaged inso- 
lation at the top of the troposphere. Katayama’s 63 percent 
(1967), and London’s 62 percent (1957) for net flux at 200mb 
are in good agreement with the 64 percent from this study. 

The net flux at the surface was 18 percent for all three 
studies. Converting these spatially averaged results to 
annual means (ly day ^) , we get a net outgoing radiation at 
P ~ 200 mb of 430 from the Kat ay ama-London means compared 
to our 450, At the surface, the net outgoing long-wave radia- 
tion is 120 (ly day ^) compared to our 131. 

In the net flux calculations, a comparison was made be- 
tween three different atmospheric situations in order to 
evaluate the contribution of cloudiness to the total outgoing 
flux from the earth-atmosphere system. One test considered 
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TABLE V(a). Net long wave flux (ly day ) distributions with latitude at reference 
levels C » 0.2, 0.6 and 1.0 with the radiative temperature cooling 
(°C day”^-) between reference levels with average total cloud cover. 
(Spring) 
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TABLE V(b). Net long wave flux (ly day ) distributions with latitude at reference 
levels a = 0.2, 0.6 and 1.0 with radiative temperature cooling (°C day 
between reference levels, with average total cloud cover. (Summer) 
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TABLE V (c) . Net long wave flux (ly day ) distributions with latitude at reference 
levels o = 0.2, 0.6 and 1.0 with radiative temperature cooling (°C day 
between reference levels, with average total cloud cover. (Fall) 
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TABLE V(d). Net long wave flux (ly day ) distributions with latitude at reference 
levels o = 0.2, 0.6 and 1.0 with radiative temperature cooling (°C day 
between reference levels, with average total cloud cover. (Winter) 
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only clear skies, another had average cloudiness and the 
third tested totally the overcast model. Results were com- 
piled in mean annual latitudinal distributions of long- 
wave net flux in Tables Vl(a-c). 

Since temperatures normally decrease with height in 

the atmosphere, a positive net flux occurred at all levels. 

Strong surface temperature inversions can lead to negative 
* 

FlO > but these conditions were not present in this study. 
Also, in an atmosphere where the effective optical depth 
decreased sharply with height, the net radiative flux in- 
creases sharply with height at the levels of the decreased 
water-vapor mass. This happens primarily because most of 
the emitted flux from the lower layer passes upwards un- 
attenuated by the dry upper layers. 

The net flux through the tropopause is larger at all 
latitudes for clear skies than for conditions of average 
or overcast cloudiness. This occurs because the net outgo- 
ing radiation comes from a lower, and therefore warmer 
layer in the case of clear skies. In the case of overcast 
skies, the radiating layer is the cloud top at a = 0.4 
Rather than a = 1.0 for clear skies. The tropopause net 
flux with average cloudiness is a weighted average between 

the overcast and clear-sky cases of Table VI. 

The largest variability of net flux magnitudes with cloud 
cover occurred at the O = 0.6 level. The minimum value of 
zero is observed during total cloudiness since this level lies 
within the cloud layer which was treated as a blackbody. Maxi- 
mum values occurred of course in conditions of clear skies. 
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TABLE VI(a). Mean annual radiative long-wave flux (ly day ) normal cloud 
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IV, SOLAR RADIATION TREATMENT 



A. EFFECTIVE SOLAR INSOLATION 

The solar radiation treatment is not as straightforward. 
Following Joseph (1966), a solar constant of 2.00 ly min ^ 
was assumed. A 4% attenuation by oxygen and ozone above 
the tropopause was assumed which gave an effective solar 
constant of 1.92 ly min The effective tropopause height 

was assumed to be at a = 0.2. 

The effective insolation at the tropopause was computed 

from 

F = S(r/r )“ 2 D cos 9 (35) 

m 

where 

S = effective solar constant 

r/r = radius vector of earth relative to the sun 
m 

0 = zenith angle of sun relative to observing point 

cos 0= daily average of cos 0 
D = fractional length of daytime 

The (r/r ) is listed as a function of the Julian date 
m 

in the Smithsonian Meteorological Tables (List, 1958). The 
following table, after Manabe and Moller (1961) was employed 
to derive the daily mean cos 0 and duration of insolation 
as a function of latitude and time of the year. Computa- 
tions of F were based upon mean mid-seasonal values of 

Spring - April 15th 
Summer - July 15th 
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TABLE VII. The seasonal and latitudinal distributions of 
duration of insolation and daily mean values of cos 0. 



Lat. Fractional Length 
of Daytime 



Apr 



Ju 1 



Oct 



Jan 



Wei ghte d Mean 
cos 

Apr 



5 


.508 


.517 


15 


.521 


.537 


25 


.533 


.562 


35 


.546 


.596 


45 


.562 


.637 


55 


.596 


.708 


65 


.629 


.837 


75 


. 750 


1 . 000 


85 


1.000 


1.000 



500 


.496 


.625 


492 


.471 


.618 


483 


.450 


.599 


471 


.421 


.558 


454 


.362 


.501 


437 


.321 


.423 


404 


.208 


.345 


329 


- 


.241 








.168 



6 



Jul 


Oct 


J an 


587 


.614 


.591 


601 


.579 


.549 


593 


.524 


.474 


567 


.458 


.393 


521 


.379 


.317 


453 


.282 


.203 


369 


.176 


.106 


311 


.071 


- 


318 
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TABLE VIII. Incoming solar radiation on a horizontal 
surface at the top of the tropopause (ly day“l). 



Latitude 


Spring 


Summer 


Fall 


Winter 


0-10 


884 


867 


844 


784 


10-20 


896 


922 


783 


692 


20-30 


889 


952 


696 


571 


30-40 


848 


965 


593 


443 


40-50 


784 


948 


473 


306 


50-60 


702 


916 


339 


174 


60-70 


604 


882 


195 


59 


70-80 


503 


888 


6 4 


0 


80-90 


468 


908 


0 


0 
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Fall 



October 15th 



Winter - January 15th 

Joseph’s method (1966) was followed in partitioning the 
effective solar insolation into a part subject to water- 
vapor absorption but not to Rayleigh scattering, and the 
remainder subject to only Rayleigh scattering. The presence 
of clouds introduced cloud reflectivities in both portions. 
Thus, the solar insolation subject to tropospheric absorp- 
tion comprises wavelengths ■> 0.9 ym and was obtained from 

F. = 0.349F (36) 

A 

In the wavelengths 0.9 ym, water vapor absorption is 
negligible and the solar insolation subject to Rayleigh 
scattering was obtained from 

F g = 0.651F (37) 

B. DISPOSITION OF FLUX SUBJECT TO SCATTERING 

In treating the fraction of the solar insolation subject 
to scattering, Joseph (1966), found that the Rayleigh 
scattering values at sea level in clear skies (after Coulson, 
1959) could be fitted by least squares to the form 

cig = 0.085 4- 0.25074 log (~ sec 9) (38) 

where 7T is the observed surface pressure. For a full-day 
the value of sec 9 in (38) was treated as sec 9, the time- 
mean sec 9, which was taken as the inverse of the mean 
cos 0 values from Table VII. 
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Another reflective parameter is the surface albedo, a . 

g 

After numerical experimentation designed to test for a 

global radiative heat balance, a constant average surface 

albedo of a = 0,14 was deduced. It should be noted that 
g 

this value of turned out to be equal to that deduced by 

Budyko (1955). A discussion of the validity of a constant 
surface albedo will be reserved for the surface balance 
section . 

Considering the possibility of multiple reflections be- 
tween earth and atmosphere, each of which diverts downward 

the fraction a of the earth surface reflectance a , the 
s g 

clear sky insolation at earth after scattering is given by 

I = F (1-a.) [1+a a + (a a ) n + ] (1-a ) (39) 

soss Sg Sg g 



or 



^SO F s ^ a g) ^ 1 -a S a g) 



(40) 



For the existence of a cloud cover, the Rayleigh re- 
flection coefficient of (38) was applied at the cloud top 
O = 0.4, rather than 0=1. At the cloud top, a unique 
cloud reflectivity fraction of 0.5 was chosen everywhere, 
This reflection is basically a Mie-type, diffuse back- 
scattering caused by water droplets in the clouds. 

Thus, at the cloud top, the entering solar-beam insola- 
tion is given by n , 

F 0 [l-a 0 (~^ sec 9 ) ] (1-R ) 

£ (41) 



SC4 






This downward scattered beam is forced into a diffuse 
angular array by the Mie scattering process. In order to 
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account for further attenuation by Rayleigh scattering when 
the ray paths are diffuse, a mean slant path angle 9^ such 
that sec 0^ = 5/3 was judged more appropriate than the actual 
sec 0 (Katayama, 1966). 

The Rayleigh scattering coefficient beneath the cloud 
top (i.e. in layer (4,10)) was found to be 

“SD - a SD,„ - “SD, - °- 25074 lo S 2 ' 5 < 42 > 

10 4 

Thus, a below the cloud top has a constant value of 0.0998. 

O U 

In this connection, it may be noted that Manabe and Stride- 
ler (1964) employed a constant value of = 0.07 S cos 0 
loss to space by Rayleigh scattering regardless of the pres- 
ence of clouds . 

Multiple reflections take place below the cloud top. 
Martin (1972) has shown that upon simplification, the inso- 
lation at earth after cloudy-sky scattering is 



r sc* F SC4 (1 -“sD H1 ' a g ) + 



F SC4 ta SD + “g (1 '“SD ) 1 (1 -V E C 



[ l-E c ( 1 -a SD ) a g ] 



(43) 



In the above accounting, certain small atmospheric re- 
flection terms occurring in the direction cloud to earth 
(or reverse) have been ignored as small compared to the air 
transmission factor (1 - ot^) . 

Finally, the weighted mean earth-surface absorption of 
the scattered solar beam becomes 



x s (l°) 



(i-CL) 1 ^ + I SC (CL) 



(44) 



where I is given in (40) and I in (43). 

j(j JU 
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C. TROPOSPHERIC ABSORPTION OF SOLAR RADIATION 

In considering the absorption of the component of 

solar insolation, the Manabe-Moller solar-absorptivity 
f unct ion 

a = 0.0946 ( u 2 2 k sec 0) 0,303 (45) 

was used. The form of Equation 45 indicates that the ab- 
sorption function a is applied over the pressure-scaled 
water vapor mass in the layer (2,2K) along the zenith 
path of slant angle 0. Restricting absorption to the ratio 
0.349F requires that the multiplicative constant 0.0946 of 
(45) be adjusted to (0.0946/0.349) = 0.271. 

The absorbed insolation in the layer (2,6) is then 
A ( 2 , 6 ) = 0.271[u(2,6)sec 6] 0 ' 303 F A (46) 

Similarly, the absorbed insolation in the layer (2,10) is 

A ( 2 , 10 ) = 0.271[u(2,10)sec e]°* 303 F A (47) 

Therefore, the absorbed insolation in the layer (6,10) is 
given by 

A (6 , 10) = A ( 2 , 1 0 ) - A ( 2 , 6 ) (48) 

Since the amount of solar beam absorbed in the layer (2,10) 
is available by (47), the solar insolation in the F^-parti- 
tion received at earth is given by F^- A(2,10). This por- 
tion is reduced by ground albedo so that the 

earth-absorbed insolation resulting from F^ is 

I. .(10)'- F. {l-0.271[u(2,10)sec 6 ] 0 * 303 } ( 1-a ) (49) 

AU A & 
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A modification is made for cloudy skies. The absorption 
in the layer (2,4) was computed according to the Manabe- 
Moller absorptivity function [Equation (46)], The solar 
beam at level 0.4 is given by 

F a ( 4) = F a - A ( 2 , 4 ) (50) 

Of the F^-energy available at level O = 0.4, the fraction 
(1-R^) is transmitted through the cloud top. Thus 

F A (4) “ [F A -A(2,4)] (1-R c ) (51) 

Within the cloud layer, the water vapor path below a = 0.4 
is subject to the mean diffuse-path augmentation factor 5/3 
(after Katayama, 1966). 

Inside the cloud layer, another variable has been intro- 
duced. The value is u^ or equivalent water vapor mass for 
droplet absorption in clouds. The values of the absorption 

coefficient a due solely to cloud particles have been 
C L iJ 

taken after theoretical computations of Korb, Michalowski 
and Moller (1956). These values are listed after Manabe 
and Strickler (1964, p. 367) in Table IX, below. 

TABLE IX. Ab s or p t i vi t i e s ( a njVi ) of cloud droplets. 

C L 13 



Latitude Band a CLD 

0-39N, low- t o-middle clouds 0.035 
40-64N, low-t o-middle clouds 0.020 
65-90N, ice-crystal clouds 0.005 



Since only a single cloud layer was used in the radia- 
tive model, the larger ab sorp t iv i t ie s associated with the 
larger water content was identified with the 0-39N latitude 
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band of Table IX. The values of a are due solely to ab- 

L* JJ 

sorption by either liquid water or the solid ice content of 
clouds. The values of a ^ are considered mean values for 
the latitude bands shown in Table IX. Moreover, the a 

U L L> 

values of Table IX are averaged for the entire solar beam, 
and not merely for the portion of the insolation. Hence, 
the value of a CLD has been multiplied by (1.0/0.349) for 
application to the F -partition. In this spectral region, 

A 

the effective cloud-water absorp tivities become 

a* = 2.8653 a CLD (52) 



or 

* r 5 , .. .,0.303 r 5 .. .. ,0.303-,- 

a = {j(u(4, 8 )+u c ) ] — [ -3 u (4,8)] }0.271 (53) 

* 

Equation (53) may be solved for as a function of a and 
u(4,8) 



u 



u ( 4 , 8 ) 



= {[ 1 + 



a ^ 0.6 ^ 0 . 303 j 3.003_ 1 



0.271u(4,8) 



(54) 



Values of u^ were computed as a function of u(4,8) for the. 
latitudes of Table IX and are shown in Appendix Table B. 
With u^ known, the following calculations are made 

A (4 , 6 ) = F A (4){0.271[| u(4,6)+ ^] 0,303 } (55) 

A(6,10)=F a ( 4) {0.271[|(u(4,10)+u c ) ] 0,303 }-A(4,6) (56) 

where F (4) is evaluated upon entering the cloud top. 

A 

Finally, the insolation reaching earth after absorption in 
cloudy skies is 

I AC (10) = F A (4) ~ A < 4 » 6 > " A( 6 , 10) (57) 
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With the surface albedo effect (l~a ) included, the earth- 

8 

absorbed part of I^^(IO) becomes 

I AC (10) = I AC (10) [ (l-O g )/(l-Cl g R c ) ] (58) 

D. NET INSOLATION ABSORBED AT EARTH 

The weighted sum for the insolation at earth after ab- 
sorption is given by 

I AC (10) = CL[I ac of (58)) + ( 1 -CL ) [ I A q of (49)] (59) 

Thus the total insolation absorbed at earth is given by 

Q l0 = I A (10) + I s (10) (60) 

The resulting values of Q^q are listed in Table X. 

E. COMPARISON OF MERIDIONALLY AVERAGED RADIATION 
QUANTITIES 

In order to compare the various annual outputs of this 
radiation model with the results of other researchers, the 
radiational quantities derived by the NPS model were time- 
averaged and space-averaged over the Northern Hemisphere. 
The resulting global mean annual quantities are compared 
with the results by Katayama (1967), Budyko (1956) and 
London (1957). The mean values are presented as a percent- 
age of the solar radiation at the top of the troposphere. 

Meridionally-averaged values were generated using the 
formulation 



R 



^2 2 

2tt / a R cos <j) d<j> 

4 >, 



^2 2 

2t r / a cos <j> d<j> 



(61) 
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TABLE X 



Latitude 

0-10 

10-20 

20-30 

30-40 

40-50 

50-60 

60-70 

70-80 

80-90 



Seasonal distribution of insolation absorbed at 
earth (Q^g) in ly day~^. 



Spring 


Summer 


Fall 


Winter 


451 


430 


424 


409 


484 


471 


404 


385 


487 


510 


374 


316 


450 


525 


315 


234 


405 


484 


241 


157 


361 


451 


166 


88 


317 


430 


95 


30 


264 


430 


31 


0 


248 


454 


0 


0 
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where 



R = annual mean radiative quantity at latitude <j> 
a = radius of earth 

<J) = center of the latitude band under consideration 

dcf> = latitude increment in degrees between successive 
bands. This increment is a constant of 10° for 
this study. 

The formulation in Equation (61) may be reduced to a 
series summation using cosine weighting factors as indicated 
in 



R 



9 



E 

n = l 



R cos <f> 
T n 



9 



E 

n = l 



cos 6 

n 



(62) 



where (j) = 5°, cj)^ = 15°, etc. 

When the results of this study are compared with those 
of Katayama, Budyko , and London in Table XI, it is found 
that they are in rather good agreement as a whole. However, 
there are some minor differences. First, the older model 
estimates of solar radiation absorbed in the atmosphere are 
generally 3. to 5 units higher than determined by this study. 
Also, while the average global albedo, A (<J>), determined in 
this study is typical of usual values quoted, namely 34%, 
the meridional distribution of A^(<J>) was more nearly uniform 
with respect to <j) than would have resulted if the means had 
been based upon hourly computations for the months involved. 
This result leads to a lower than normal value of solar 
radiation absorbed at the surface in low latitudes and a 
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TABLE XI. Comparison of meridionally averaged values of 
radiation quantities. The values are presented as per- 
centages of the global mean insolation at the top of the 
troposphere. 

Katayama Budyko London Plante 



Solar radiation at 
top of troposphere 


100 


ioo 


100 


100 


Solar radiation absorbed 
in atmosphere 


18 


16 


14 


13 


Global albedo 


37 


40 


35 


34 


Solar radiation absorbed 
at surface 


45 


44 


47 


52 


Outgoing long -wave 
radiation at 200 mb 


63 


60 


62 


64 


Outgoing long -wave 
radiation at surface 


18 


17 


18 


18 


Sensible and latent 

heat flux (see Section IX) 


29 


27 


29 


31 
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larger than typical value at the surface in high latitudes. 
The reason for this is basically threefold: First, atmos- 

pheric absorption by dust was not considered in the NPS 
model, since it is essentially a predictive model. However, 
London (1957) applied the results of Houghton (1954) who 
used atmospheric transmissivity due to dust as given by 

t d - ( 0 . 9 5 ) m (63) 

where m - sec 0, represents the air mass traversed by the 
solar beam. London’s depletion by atmospheric dust would 
then be given by 

D « <1 - t d ) (64) 

of which 25% is considered subject to absorption. Even 
so, the contribution of dust absorption is small and can 
generally be considered of the order of 1%, which is 
approximately the difference between the NPS absorption 
figure and that of London’s. 

Another contribution to absorption not previously con- 
sidered here is secondary absorption of the reflected por- 
tion of the solar beam. Since the reflected portion of the 
solar beam above clouds is subject to the same percentage 
absorption as the incoming beam, the error in absorption 
may be as high as 1% of the incoming solar radiation. Neither 
London’s study nor this study considered secondary absorp- 
tion of the reflected beam and consequently the resulting 
atmospheric absorptions are closely comparable. Katayama 
considered secondary absorption following a first reflection 
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of the solar beam by a cloud or by ground or both. As a 
result Katayama’s solar absorption in air was the largest 
of the four cases considered. 

The third consideration is Katayama’ s cloud cover. In 
the Katayama radiative study, a larger climatological cloud 
cover was considered than in the present study, or that of 
London or Budyko . The cloud amounts were generally charac- 
teristic of coverages over the oceans (after McDonald, 1938) 
rather than a combined ocean-continent regime. 

It is not exactly clear how Katayama allowed McDonald’s 
overlapping cloud-layers to be combined, although he con- 
strained the total cloud-coverage to agree with that of Sir 
Napier Shaw’s (1936). The significant feature of the 
greater cloud coverages is that they afford not only greater 
solar absorptivity but also a greater global albedo as 
postulated by Katayama’s empirical formulas 



where h is the cloud thickness and b^, b^y are 

characteristic values for a given cloud type. Thus Kata- 
yama finds a greater cloud absorptivity than was found by 
London or Plante (Table XI) as well as a larger global 
albedo because of the larger cloud coverages. 

F. GLOBAL ALBEDO CONSIDERATIONS 

In the process of deriving the transmission of the two 
solar components F A and F , the following dispositional 
terms were discussed: 



A c = [ 1-exp ( -b 1 h) ] 

R c = R oo [ 1-exp ( -b 2 h)] 



(65) 



( 66 ) 
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(i) the combined insolation rate at the tropopause 
denoted by F of Equation (35), 

(ii) the absorption rate of insolation in the atmos- 
phere 

(iii) the absorption rate of insolation, Q^q> at earth 

With these terms, the reflected insolational rate from 

planet earth, as a function of latitude and season may be 
computed. This rate may be denoted by A^F, where is the 

global albedo. The term A^F must be identified as the 
difference of 

F - A ( 2 , 10 ) - Q 1q = A p F (67) 

It should be noted that three significant reflective 

parameters have affected the value of the left side of (67) 

namely the sky albedo, the cloud albedo (R^ = 0.5), and the 

earth-surface albedo, a . For simplicity, the latter two 

& 

were held constant at reasonably well accepted figures in- 
dependent of latitude. While direct computations of the 
right side of (67) were made at each step of computing (ii) 
and (iii), the resulting global albedo values must also be 
consistent with the values inferred from the left side of 
(67) . 

A note of caution in regard to comparing global albedo 
as shown in Table XI. A correct mean global albedo may 
distort the atmospheric heat balance if not distributed 
correctly with respect to latitude. For example, if a 
meridional bias exists in the computed global albedo, A^ , 
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with too large a value at low latitudes and too small a 
value at high latitudes, the result will influence the pole 
ward heat transport required for balance, after independent 
computation of latent and sensible heat transfers across 
the earth-air interface (see Section IX) • This type of 
bias in global albedo appears to have developed in this 
study presumably by oversimplifying the ray-paths in clouds 
Choosing a single mean path through clouds for a given re- 
presentative day, and time of year seems to have overesti- 
mated the global albedo in this study in lower latitudes 
while underestimating it at high latitudes. 
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V. RADIATIVE TEMPERATURE CHANGES 



Long-wave radiation, with a normal temperature lapse, 

will lead to a cooling of the atmosphere as seen as seasonal 

distributions in Tables v(a-d). In certain locations and 

seasons, some of the atmospheric layers considered were 

subject to temperature inversions which led to local heating 

rather than the usual cooling. Another heating situation 

examined was at an overcast cloud base. The terminology, 

’cloud base warming’ or ’green house effect’ may be used 

to describe this phenomenon. Since the net flux within a 

cloud layer is zero, a negative net flux increment existed 

in the cloud to sub-cloud layer (8,10). With the sign 

convention adopted here, this represents a base-warming 

effect. In the cloud top region, there is a substantial 

* * 

net flux increment F 2 - F^ which leads to considerable 
cloud-top cooling. Thus, in infrared radiative temperature 
changes the role of clouds is clearly to increase the cooling 
rate at the cloud tops and to decrease the cooling rate 
at the cloud bases. This process can be considered as a 
decrease of the atmospheric stability within the cloud layer. 

The mean annual infrared temperature changes in the 
troposphere are plotted in Figure 7 as a function of latitude. 
The dependency of these cooling rates upon the climato- 
logical cloud cover at these same latitudes is evident by 
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inspecting Figure 8. A linear regression was performed re- 
lating the mean long-wave temperature changes of each tropo- 
sphere layer against the total cloud cover used in the 
computation at latitude (p . Each latitude band for each 
season was considered separately, giving a total of 36 cases. 
For the layer (6,10) the correlation coefficient was -0.95 
while for the layer (2,6) the correlation was 4-0.99. Graphs 
showing the comparison of temperature changes found with 
clear and with average cloudy skies for the layer (2,6) and 
layer (6,10) are shown in Figures 9 and 10, and clearly 
demonstrate the phenomena of cloud-top cooling and cloud- 
base warming. 

In conditions of normal cloudiness, the infrared heat 
losses are partially offset by atmospheric absorption of 
solar radiation. The resultant temperature change in the 
layer (6,10) is 



st ~(F *- F in* ) + A ( 6 > 10 ) 

9T. /£ nr >\ _ o 10 

3 1 ^ ' “ 0 . Alt (1000/ g) c 



( 68 ) 



The temperature change in the layer (2,6) is given by 
ST -(?,*- F *) + A ( 2 , 6 ) 

ft (6 » 10) = 0 . 4tt (1000/g) " Cp (69) 

In most cases, this absorption of solar insolation is smal- 
ler than the cooling that results from long wave radiation 
and acts only to modify the pattern of radiative heat loss, 
as shown in Figure 11. The dashed curves represent the 
annual mean rate of temperature change due to long wave 
radiative cooling for the layers (6,10) and (2,6). The 
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Figure 11. Comparison of mean annual atmospheric temperature chang 
for infrared cooling and net radiation for layers (6,10) and (2,6). 



solid curves represent the result of net radiative cooling 
(combined effects of long and short wave radiation) in the 
same layers. Atmospheric absorption clearly only modified 
the magnitude of radiative cooling while the shape of the 
curves was unchanged. By comparison of Figures 8 and 11, 
high correlations between net radiative temperature changes 
and total cloud cover may once again be noted. 
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VI. RADIATIVE BALANCE AT THE TROPOPAUSE 



A radiative budget was calculated at the tropopause and 
compared to the results of Houghton (1954). Houghton made 
a comparison of the net incoming solar radiation at the 
tropopause with the outgoing long-wave flux. He determined 
that insofar as the earth-atmosphere combined system is con 
cerned, a column at about 35°N latitude experiences a mean 
radiative energy balance. Corresponding latitude belts of 
the system located equatorward of latitude 35N have a mean 
annual radiative energy surplus. Those latitude bands pole 
ward of 35°N latitude experience a mean annual radiative 
energy deficit. These mean annual radiative surpluses and 
deficits impose requirements for horizontal heat transfer 
rates across latitude walls in both air and ocean. 

The net insolation at the tropopause (Q q ) is given by 



a 4% reduction due to attenuation in the stratosphere. Q 
is the reflected insolation returned to space. The ratio 




(70) 



where Q^. is the incident insolation at the tropopause after 



of Q / £ 1 . 0 4 Q ] is the measure of the earth atmosphere 
R I 



planetary albedo. Q R is given by 

= ‘Ho 



(71) 



where Q is the short wave radiation absorbed in the atmos 
A 

phere and is the total short wave radiation absorbed 

at earth. 
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The resultant Q R is a function of the Rayleigh back 
scattering, cloud reflectivity and surface albedo. The 
factors determining Rayleigh scattering were discussed 
in Section IV. 

Treatment of cloud reflectivity was limited here by 
the use of a single cloud layer. A global average of 
about one-fourth of the incident solar radiation is re- 
flected to space by clouds. In this model a single cloud 
albedo R^ = 0.5 was used. Sellers (1965) employed a range 
of albedos for differing cloud types. A similar treatment 
using various cloud types with various albedos was adopted 
by London. On a mean annual basis, 0.5 appears to be a 
good approximation for the hemispheric cloud reflectivity, 
but on an individual gridpoint basis the cloud albedos 
adopted in this study may be an oversimplification. Future 
testing in this area is envisioned for the NPS predictive 
model . 

A single constant surface albedo of 0.14 was chosen 
after experimentation with the tropospheric energy balance. 
Budyko (1955) estimated the mean earth surface albedo as 
0.14 while Katayama (1967) found a value of 0.13. But the 
albedo varies with the type of surface. The short-wave 
albedo is highest over fresh snow ( 90%) and is lowest over 
water and forests (approximately 5%) . Since there are 
seasonal as well as annual variations in the extent of 
vegetation and snow cover, an attempt at strict evaluation 
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of ground albedo for these variations is at best difficult. 
Thus, for simplicity, a single mean ground albedo was em- 
ployed. 

The net outgoing long-wave radiation at the tropopause 
* 

was given by F 2 since the 0.2 sigma level was adopted as 
the mean height of the tropopause. 

Therefore, the net radiative flux across the tropopause 
is given by 



* 



Q N - Q c - F 2 (72) 

Results for seasonal variations in the net flux across the 
tropopause are listed in Tables Xll(a-d). The mean annual 
radiative input by latitude bands is listed in Table Xll(e) 
along with the cosine-latitude weighting factors. The re- 
sulting mean Northern Hemisphere radiative budget is zero 
as anticipated, since it is assumed that the hemisphere is 
in annual radiative balance. 

During the summer, the tropopause has a net radiative 
surplus throughout all latitude bands. As seen in Table 
XII(b), the incident insolation is largest during this sea- 
son and far exceeds the small increase in net long wave 
radiation due to warmer summer temperatures. The opposite 
case is true for winter. Low insolation rates produced a 
negative radiative budget at all latitude bands except the 
0-10° band. This particular latitude band has a positive 
radiative net surplus during the entire year. Spring and 
fall show transitions between the extremes of summer and 
wint er . 
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TABLE XII(a). Net radiative budget (ly day ) at the tropopause with average cloud 

cover. (Spring) 
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TABLE XII(b). Radiative budget (ly day ) at the tropopause with average cloud 

cover. (Summer) 
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TABLE XII(c). Radiative budget (ly day ) at the tropopause with average cloud 

cover. (Fall) 
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TABLE XII(d). Radiative budget (ly day ) at the tropopause with average cloud 

cover. (Winter) 
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TABLE XII(e). Mean annual radiative budget at the tropopause. 
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Effects of variations in cloud cover on the computed 
net radiative input at the tropopause can be seen in Table 
XIII which lists spatial mean averages of the radiative 
quantities across the tropopause. It is seen that clouds 
exert an overall cooling effect in this model. 



TABLE XIII. Spatial mean annual net radiative input at 
the tropopause (Northern Hemisphere). 





CL = 0 


CL = 1.0 


Average CL 


Q I 


694 


694 


694 


q r 


130 


360 


244 


Q o 
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334 
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F 2 
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373 
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Net Tropopause 
Radiative Input 
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-39 


0 
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VII. GLOBAL ALBEDO 



The earth’s global albedo results from the upward re- 
flection by the ground, by clouds and by scattering by 
atmospheric molecules of solar radiation. Variations in 
computed planetary albedo are strongly dependent on total 
cloud cover. The distribution with latitude and season 
of global albedo is listed in Table XIV. There is a pro- 
nounced minimum in the subtropics, which is a region of 
minimum cloudiness. The minimum shifts equatorvard during 
the winter and poleward during the summer, following move- 
ment of the ITCZ. The increased albedo in high latitudes 
during all seasons results from an increase in cloudiness 
with latitude up to about 60 to 70N. 

Total global albedo when averaged over all latitudes 
and seasons is about 34 percent. The total global albedo 
shows little variation with season and agrees closely with 
London’s results. Figure 12 shows a comparison between 
mean annual planetary albedo with latitude computed by 
London and the results of this study. Even though the 
spatially averaged annual global albedo of 34% agrees close- 
ly with London, there is disagreement in its distribution 
with latitude, a fact referred to in Section IV. 

In tropical latitudes, London’s global albedo was less 
than that computed in this study, while in the northern 
latitudes, his values exceeded those computed here. Since 
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Comparison of mean annual global albedo distributions with 
latitude . 



TABLE XIV. 



Seasonal distribution of planetary albedo with 
latitude (in percent of Qj)> 
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35 
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identical total cloud cover amounts were used in both 

studies, the treatment of atmospheric and cloud albedos 

seem to be the major cause for this difference. The 

present study used a constant cloud albedo of R = 0.5 

c 

while London used varying cloud albedos associated with 
constituent cloud types after Haurwitz (1948). London 
thus used a multi-layer cloud model while this was not 
done here. Furthermore, the device of replacing the 
diffuse path length in clouds by an equivalent path 5/3 
longer seems to give too short a path subject to cloud- 
reflecting processes at high latitudes and too long a 
path at low latitudes. This is important if cloud-reflec- 
tance is to be treated by Mie-scat tering theory. Varia- 
tions of surface reflectance with solar zenith angle also 
plays a secondary role in explaining the variations of 

A with latitude . 

P 
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VIII* ATMOSPHERIC BUDGET 



A simplified annual atmospheric heat budget for the 
Northern Hemisphere, along with seasonal distributions of 
radiative surpluses and deficits was constructed. The 
stratosphere was assumed to be in radiative equilibrium. 

The atmosphere radiation balance is shown schematically 
below : 



q r 

A i 

Q 10 q e q s 

I t - t 





The net cooling contribution due to long wave radiation is 
given by 

F N ■ F 2 - F 1*0 (73) 



The effective insolation at the top of the troposphere 
is the total incident solar radiation minus the reflected 
insolation to space. The total atmospheric absorption is 
therefore 



” q r- ) °ho 



(74) 



where Q is the short wave radiation absorbed at earth 

using a mean surface albedo, a = .14. A final heating 

8 

term consists of the latent heat of evaporation and sen- 
sible heat transfer rates across the earth-air interface 



and denoted by Q„ : 

H 



+ Q ! 



(75) 
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The parameterization of Q £ and Q s is detailed in Section IX. 
For each latitude band, the net atmospheric heating rate by 
the three mechanisms (insolation absorption, long-wave 
absorption, interface convection) is given by: 



{ VA 



= Q 



H 



+ - 



N 



(76) 



The resulting atmospheric heating (cooling) rates in ly day"* 
are summarized in Tables XV(a-d) for each latitude band and 
each of the four seasons, while Table XV(e) gives the annual 
mean, derived by averaging the four seasonal heating rates, 
as a function of latitude. 

It is to be noted that Q has been termed the '.'balance 
requirement" by Davis (1963), and represents the excess 
amount to be disposed of from a latitude belt by poleward 
flux-divergence of heat. Table XV(e) shows that the global 
annual average of this parameter (as calculated in this 
study) is zero. 

Warming of the atmosphere by the absorption of solar 
radiation is less than the long-wave cooling rate. The net 
result is that the atmosphere is cooling radiatively at all 
latitudes at a rate of about 1.1°C day The compensating 

energy requirement for possible balance has been indicated 
by Equation (7 6) as arising by sensible and latent heat trans- 
fers across the earth surface. 

With the addition of sensible and latent heat into each 
latitude band, Table XI(e) shows that latitude bands equa- 
torward of AON have a total energy surplus while latitude 
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bands poleward of AON exhibit an energy deficit. While the 
mechanisms of horizontal transport are not specified in 
this study j the resultant glob ally-we igh t ed values of 
of Table XV(e) indicate a mean annual total heat balance 
in the atmosphere. The poleward redistribution required to 
effect this total heat balance is evidently brought about 
by general circulation processes, which were not considered 
in this study. For our purposes it suffices to show that 
the energy-generating aspects of the model proposed here is 
subject to an atmospheric balance. 
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TABLE XV(a) . Atmospheric heating rates (ly day ), based upon climatological cloud 
cover. (Spring) 
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TABLE XV(b). Atmospheric heating rates (ly day ), based upon climatological cloud 
cover. (Summer) 
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TABLE XV(c). Atmospheric heating rates (ly day ), based upon climatological cloud 
cover. (Fall) 
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TABLE XV(d). Atmospheric heating rates (ly day ), based upon climatological cloud 
cover. (Winter) 
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TABLE XV(e). Mean annual atmospheric heating rates (ly day 
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IX. SURFACE HEAT TRANSFERS 



The earth's surface heat balance was computed using 
evaporation and sensible heat transfers over the oceans 
as surface heat sinks. Heat transfer due to evaporation 
and sensible heat are more easily computed over ocean 
surfaces using climatological data available in the Paci- 
fic Ocean representative of the latitude belt as a whole. 
In effect, London's atmospheric soundings were used in 
association with data for ocean surface temperatures and 
surface winds taken from U. S. Weather Bureau's Atlas of 
Climatic Charts of the Oceans (1938). Continuity of 
sounding data with the oceanic reports was maintained by 
maintaining the same temperature lapse rate in the layer 
zero to 1 km. 

The surface radiative budget is shown schematically 
below : 





Q 



s 



1 1 



and may be stated symbolically by the equation: 

Q io - (F io* + q e + V ■ «vo 



(77) 



Here Q vq is the required oceanic heat transport divergence 
The surface layer evaporation heat transfer rate was calcu 
lated after Langlois and Kwok (1969) using: 



Q E ~ L P 10 C D V S 



sg q x^ 



(78) 



88 



with the constants: 



P 



L 

10 



590 cal gram 
-3 

1*2 x 10 gram cm 



Also V is the climatological ocean surface wind speed in- 
cremented by 2.2 mps in order to account for the effect 
of gustiness upon evaporation [Langlois and Kwok (1969)]. 

A variable drag coefficient proposed by Deacon and Webb 

(1962) for the case of neutral stability was employed: 

C = (1.0 + 0.07 V ) x lo" 3 (79) 

JJ 



At the ocean surface, the atmosphere was assumed to 
be saturated with respect to a plane water surface. The 
surface specific humidity, q^g, was computed using the 
Clausius- Cl apeyron equation: 



. 622 e (AE-BE/TG) 
q sg n-e (AE-BE/TG) 



(80) 



The specific humidity at the top of the constant flux sur- 
face layer, q , is given following the parameterization of 

X 

the planetary boundary layer by Langlois and Kwok, by the 
formula : 



with: 



q x = 



kq 9 + ( ' Z 9~ Z lo' ) C D V S q sg 



k + ( V Z 10> C D V S 



k 



1 + 



* 

k 



a* ( 6 9 " 9 s ) 




(81) 



(82) 



In (82), k and a are turbulent transfer coefficients 
whose values have been documented by Kaitala (1972). 
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Ice free oceans exist equatorward of 60N for all seasons 
Therefore, air and sea temperature and wind data were re- 
quired for only the first six latitude bands. Since Lon- 
don’s atmospheric data was distributed at one-kilometer 
intervals, (Zg - Z^^) was taken as 1 km; also Tg and qg are 
temperature and specific humidity at the one-km height. The 
first kilometer temperature lapse rate was computed using 
London’s soundings. Air temperature at the ocean surface 
was assumed equal to the ocean surface temperature. The 
lapse rate was then applied to yield a new T n for use in 

y 

(82). The new specific humidity q at 1 km was generated 

y 

using the relative humidity from London’s data and the 
Clausius Clapeyron equation. 

Sensible heat transfer rates equatorward of 60N were com 
puted using Bowen ratio (r) values from Sellers (1965) and 
the formulation: 

Q s - r Q e (83) 

Poleward of 60N, combined evaporation and sensible heating 
rates were taken from an assessment of many thousands of 
ship (and also of ice island) reports made by Vowinckel and 
Taylor (1964) . 

Tables XVl(a-d) contain the seasonal surface radiative 
balances, Q v0 > which is the oceanic heat transfer require- 
ment by poleward flux— divergence , Seasonal variations can 
be compared against the annual mean values listed in Table 
XVI (e) . 
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The physical meaning of a negative Q value for a 
latitude belt in the tropics (with positive values in 
northerly belts) is that a flux-convergence (divergence) 
of sensible heat due to ocean currents is needed to 
maintain a long-term balance. 
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TABLE XVI (a). Surface heating rates (ly day ). (Spring) 
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TABLE XVI(b). Surface heating rates (ly day ). (Summer) 
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TABLE XVI(c). Surface heating rates (ly day ). (Fall) 
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TABLE XVI (d) . Surface heating rates (ly day )• (Winter) 
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TABLE XVI(e). Mean annual surface heating budget estimated by the NPS model 

(ly day -1 ) . 
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X. CONCLUSION 



With the use of mean seasonal atmospheric data, a pro- 
posed radiation-heating model for the Naval Postgraduate 
School primitive-equation predictive system was evaluated. 
The model proved readily adaptable to computerization. 

The time and space averaged radiational quantities 
showed rather good agreement as a whole when compared with 
the results of other researchers. Minor difference in the 
comparisons can generally be ascribed to the treatment of 
global albedo, and the assumption of but a single cloud 
layer. While the average global albedo determined in this 
study is typical of usual values quoted, namely 34%, the 
meridional distribution of global albedo was more nearly 
uniform with respect to latitude than would have resulted 
if the mean value had been based upon hourly computations 
for the months involved. This result led to a lower than 
normal value of solar radiation absorbed within the earth’s 
surface in low latitudes and a larger than typical value 
absorbed in high latitudes. 

Of the three heat budgets discussed, results compatible 
with Houghton (1954) and Davis (1963) respectively, were 
demonstrated for the earth-atmosphere system and for the 
magnitude and direction of the required total heat flow. 
Similar tests were also applied to the ocean surface 
where the distribution of surface heat-flux divergence was 
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opposite to that expected (e.g. Budyko ( 1955 )), the results 
of the surface heat budget were somewhat biased by the 
averaging involved in the derivation of the global albedo, 
as explained previously. 

Finally, while the radiation package has proven general- 
ly adaptable to a general circulation predictive model, fur- 
ther refinements are anticipated in the treatment of cloud 
cover, A two-layer cloud model for the treatment of middle 
and low-level clouds has been suggested. These two cloud 
types would have different short-wave flux reflection and 
absorption characteristics and would lead to a more accurate 
determination of global albedo. 
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APPENDIX A 



.Evaluation of the final temperature-dependent flux integral 
e (u , T^) dB for a variety of final temperature values (T^) 
and final scaled water-vapor values (ly day~^) 



Scaled 



Water 

Vapor 


190° 


195° 


200° 


205° 


210° 


215° 


220° 


0.0001 


26.98 


28.53 


30.12 


31.76 


33.45 


35.19 


36.97 


0.0005 


42.62 


45.36 


48.19 


51.13 


54.16 


57.30 


60.55 


0.001 


50.84 


54.25 


57.80 


61.48 


65.30 


69.26 


73.36 


0.002 


59.89 


64.09 


68.46 


73.01 


77.75 


82.66 


87.78 


0.004 


69.65 


74.74 


80.05 


85.60 


91.38 


97.41 


103.68 


0.006 


75.63 


81.29 


87.20 


93.39 


99.85 


106.59 


113.61 


0.008 


79.97 


86.05 


92.41 


99.08 


106.04 


113.32 


120.91 


0.01 


83.37 


89.79 


96.52 


103.57 


110.94 


118.64 


126.69 


0.015 


89.63 


96.68 


104.09 


111.86 


120.00 


128.53 


137.44 


0.02 


94.09 


101.61 


109.52 


117.83 


126.54 


135.67 


145.22 


0.03 


100.38 


108.57 


117.21 


126.30 


135.83 


145 . 84 


156.34 


0.04 


104.80 


113.49 


122.65 


132.30 


142.45 


153.11 


164.29 


0.06 


110.94 


120.33 


130.25 


140.71 


151.73 


163.33 


175.51 


0.08 


115.20 


125.09 


135.56 


146.60 


158.25 


170.52 


183.42 


0.1 


118.42 


128.71 


139.60 


151.10 


163.24 


176.04 


189.51 


0.15 


124.08 


135.08 


146.73 


159.07 


172.11 


185.87 


200.37 


0.2 


127.90 


139.40 


151.59 


164.52 


178.19 


192.64 


207.89 


0.3 


132.96 


145.15 


158.10 


171.84 


186.40 


201.81 


218.10 


0.4 


136.30 


148.96 


162.43 


176.74 


191.92 


208.00 


225.01 


0.5 


138.72 


151.74 


165.61 


180.35 


196.00 


212.60 


230.16 


0.6 


140.59 


153.89 


168.07 


183.16 


199.19 


216.20 


234.22 


0.8 


143.31 


157.05 


171.72 


187.34 


203.95 


221.59 


240.30 


1.0 


145.22 


159.29 


174.32 


190.33 


297.38 


225.50 


244.73 


1.5 


148.23 


162.86 


178.50 


195.20 


213.00 


231.96 


252.10 


2.0 


149.99 


164.97 


181.02 


198.17 


216.48 


235.98 


256.73 


2.5 


151.12 


166.36 


182.70 


200.18 


218.85 


238.76 


259.97 


3.0 


151.89 


167.33' 


183.90 


201.63 


220.59 


240.81 


262.36 


4.0 


152.82 


168.55 


185.44 


203.55 


222.92 


243.61 


265.68 


5.0 


153.29 


169.23 


186.35 


204.72 


224.38 


245.41 


267.86 


6.0 


153.53 


169.61 


186.90 


205.46 


225.35 


246.63 


269.40 


8.0 


153.59 


169.88 


187.41 


206.24 


226.45 


248.09 


271.23 


10.0 


153.40 


169.81 


187.50 


206.51 


226.94 


248.82 


272.25 


15.0 


152.49 


169.08 


186.98 


206.26 


227.00 


249.26 


273.11 


25.0 


150.35 


167.04 


185.09 


204.58 


225.57 


248.14 


272.38 


30.0 


149.32 


166.02 


184.09 


203.61 


224.66 


247.30 


271.63 
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APPENDIX B 



Water-vapor equivalent u for droplet absorption in relation 
to the pressure-scaled water-vapor absorber mass u(4,8) 
(in the sounding layer (4,8)). Case u (1) corresponds to 
low latitude water-drop absorptivity; C u (2) to mid-latitude 
absorptivity; u^(3) to high latitude absorptivity. 



u (4 , 8) 


u (1) 
c 


u (2) 
c 


u (3) 
c 


0 . 01 


0 . 14170 


0 . 05118 


0.00740 


0 . 02 


0 . 18936 


0 .97292 


0.01155 


0.05 


0.28886 


0 .12003 


0 .02103 


0.10 


0 . 40852 


0 .17851 


0 .03329 


0.15 


0.50529 


0 . 22673 


0.04365 


0 . 20 


0.56995 


0.26939 


0.05294 


0.30 


0.73785 


0.34473 


0 .96955 


0.50 


0 .98617 


0.47287 


0.09826 


0 . 80 


1.29740 


0 . 63548 


0 . 13520 


1.00 


1.48119 


0 .73226 


0.15739 


1.20 


1 .65219 


0.82271 


0 .17822 


1.40 


1.81232 


0 .90822 


0.19800 


1 . 60 


1.96648 


0.98974 


0 . 21691 


1.80 


2.11306 


1.06793 


0 .23510 


2 . 00 


2 .25403 


1 . 14331 


0 . 25267 


2 . 20 


2 . 39016 


1.21623 


0 . 26970 


2 .40 


2 .52206 


1.28699 


0 .28626 


2 .60 


2 .65022 


1.35584 


0.30240 


2.80 


2.77499 


1 .42297 


0 .31816 


3.00 


2 . 89675 


1.48854 


0.33357 
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APPENDIX C 



The distribution of relative humidity and total optical 
mass (to the tropopause) u(gm cm -2 ) (from London, 1957). 
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APPENDIX D 



The distribution of average temperature and pressure (from 
London , 1957 ) . 
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T(*C) p(«b) t<*C) p(ab) 7(*C) p(xb) T(*C) p(ab) ' TCC) p(ab) T(*C) p(sb) T(*C) ?(»b) T(*C) p<ab) T(*t) p<cb3 

Q.io-- io-2o*s zzixi i *q-5q*s 50-60*7? 6o-?o*y 

Sfc 27.3 ion 27.9 1011 26.9 1012 23.2 1C1* 13.3 1013 13*3 1011 9.6 1011 3.* 101* C.* 1C13 
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1C11 27.7 1012 25.5 1016 20.1 1013 10.7 1X7 3.9 1013 - 4.4 1011 - 7.4 1013 -14.3 1015 
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p(*b) T ( *C) p(»*) 
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